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C3 i 1. Introduction 



We devote this introductory section to a historical account of some relevant contri- 
butions to the subject of operator zeta functions and their applications in theoretical 
physics. We start recalling the definition of the standard zeta function: the Riemann 
zeta function, £(s), which is a function of a complex variable, s. To define it, one 
considers the infinite series 

n=l 

which is absolutely convergent for all complex values of s such that Re s > 1, and 
then defines £(s) as the analytic continuation, to the whole complex s— plane, of the 
function given for Re s > 1 by the sum of the preceding series. Actually, Leonhard 
Euler had already considered the above series in 1740, but only for positive integer 
values of s and some years later Chebyshev had extended the definition to Re s > 1. 
Riemann formulated his famous hypothesis on the non-trivial zeros of the zeta 
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function in 185<ffl It turns out that this is the only one in the famous list of twenty- 
three problems discussed in the address given in Paris by David Hilbert (on August 
9, 1900), which has gone into the new list of seven Millennium Prize Problems, 
established by the Clay Mathematics Institute of Cambridge, Massachusetts (USA), 
and which were announced at a meeting in Paris (held on May 24, 2000) at the 
College de France. 

In 1916, in their seminal paper "Contributions to the Theory of the Riemann 
Zeta-Function and the Theory of the Distribution of Primes ''^ , Godfrey H. Hardy 
and John E. Littlewood did much of the earlier work concerning possible applica- 
tions of the zeta function as a regularization procedure, by establishing the conver- 
gence and equivalence of series regularized with the heat kernel and zeta function 
regularization methods. Also very important in this respect was the appearance 
of Hardy's book entitled Divergent Series^. We should also note that Srinivasa I. 
Ramanujan had already found, working in isolation, the functional equation of the 
zeta function, independently of all this development, as Hardy could later certify. 

Torsten Carleman, in 1935, in his work in French Proprietes asymptotiques des 
fonctions fondamentales des membranes vibrante^, obtained the zeta function en- 
coding the eigenvalues of the Laplacian of a compact Riemannian manifold for the 
case of a compact region of the plane. This was an important first step towards 
extending the concept of zeta function as associated to the spectrum of a differen- 
tial operator, which is actually the situation at issue here. And, as has been much 
more widely recognized in the literature, a decade and a half later Subbaramiah Mi- 
nakshisundaram and Ake Pleijel, in their 1949 paper Some properties of the eigen- 
functions of the Laplace- operator on Riemannian manifold^ extended Carleman's 
results explicitly showing that, if A is the Laplacian of a compact Riemannian man- 
ifold, then the corresponding zeta function, Ca(s), converges and has an analytic 
continuation as a meromorphic function to all complex numbers, what is actually a 
very remarkable result. 

Another milestone in this development was Robert Seeley's seminal work, pub- 
lished in 1967, Complex powers of an elliptic operato^. Seeley fully extended in 
this paper the above treatment to general elliptic pseudo-differential operators on 
compact Riemannian manifolds. He proved that, for all such operators, one can 
rigorously define a determinant using zeta function regularization. In 1971, Daniel 
B. Ray and Isadore M. Singer^ used Seeley's theory in their famous paper entitled 
R-torsion and the Laplacian on Riemannian manifolds to define the determinant of 
a positive self-adjoint operator, A. Such operator is, in their explicit applications, 
the Laplacian of a Riemannian manifold; denoting its eigenvalues by a\, 02, then 
its zeta function is formally given by the trace 



( A (s) = Tr (A- s ). 



(2) 
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The method defines also the (possibly divergent) infinite product as 

oo 

I] a n = exp[-a'(0)]. (3) 

At this point we arrive, in our chronological survey, to the very important contri- 
bution of Stuart Dowker and Raymond Critchley. In their seminal work, published in 
1976, Effective Lagrangian and energy-momentum tensor in de Sitter space^ , these 
authors went definitely further in the application of the above procedures to physics: 
they actually proposed, for the first time, a fully-fledged zeta function regularization 
method for quantum physical systems. This paper has got high recognition, having 
gathered over 600 citations to present date. For the sake of completeness let us here 
reproduce in full its abstract: 

The effective Lagrangian and vacuum energy-momentum tensor < T^ v > due to 
a scalar field in a de Sitter space background are calculated using the dimensional- 
regularization method. For generality the scalar field equation is chosen in the form 
(□ 2 + £i? + m )<p = 0. If £ = 1/6 and m = 0, the renormalized < T^ v > equals 
<7 A1!y (960'7r 2 a 4 )~ 1 , where a is the radius of de Sitter space. More formally, a general 
zeta-function method is developed. It yields the renormalized effective Lagrangian 
as the derivative of the zeta function on the curved space. This method is shown 
to be virtually identical to a method of dimensional regularization applicable to any 
Riemann space. 

One thing specialists often point out is that, in spite of the fact that, elaborating 
from the methods developed in this paper, it is true that a well defined and clear 
regularization prescription for a general case can be easily obtained, the authors ac- 
tually described the method very briefly in this work, the uses and wide possibilities 
of the procedure not having been fully exploited or even anticipated there. 

This is maybe the main reason why Stephen Hawking's 1977 extremely influ- 
ential paper (it has got over 1100 citations up to date) entitled Zeta function reg- 
ularization of path integrals in curved spacetim^ is considered by many to be the 
actual seminal reference where the zeta function regularization method was defined, 
with all its computational power and possible physical applications, which were very 
clearly identified there. Needless to say, the title of the paper is absolutely explicit. 
Again, let us reproduce, for comparison, its abstract: 

This paper describes a technique for regularizing quadratic path integrals on a 
curved background spacetime. One forms a generalized zeta function from the eigen- 
values of the differential operator that appears in the action integral. The zeta func- 
tion is a meromorphic function and its gradient at the origin is defined to be the 
determinant of the operator. This technique agrees with dimensional regularization 
where one generalises to n dimensions by adding extra flat dimensions. The gener- 
alized zeta function can be expressed as a Mellin transform of the kernel of the heat 
equation which describes diffusion over the four dimensional spacetime manifold in 
a fifth dimension of parameter time. Using the asymptotic expansion for the heat 
kernel, one can deduce the behaviour of the path integral under scale transformations 
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of the background metric. This suggests that there may be a natural cut off in the 
integral over all black hole background metrics. By functionally differentiating the 
path integral one obtains an energy momentum tensor which is finite even on the 
horizon of a black hole. This electromagnetic tensor has an anomalous trace. 

In my view, after investigating the case in some detail, it is fair to conclude 
that the priority of Dowker and Critchley in this matter has been sufficiently well 
established in the literature I have consulted (with some really incredible exceptions, 
however, as the running Wikipedia article on "Zeta function regularization" , where 
not the least reference to Dowker and Critchley is done!). Further to this, considering 
the number of citations collected by each one of the two papers, taking then into 
account the enormous mediatic impact of S.W. Hawking to modern physics (and 
well beyond it), and also the careful analysis of both the abstracts and the whole 
papers themselves, the ratio of citations to both works seems fair enough. But this 
is just to be taken as my personal opinion, of course. 

To continue this account further would require a very hard work and would end 
in a very long report, at least book size, what is not the purpose here. Let us finish 
this short report here, at the point when, as already mentioned, the zeta function 
regularization method is considered to have been clearly defined and its usefulness 
for physics undoubtedly established. A large number of very interesting research 
articles in several directions have been published on these matters during the last 
35 years. I will just mention a few references |10j . which by no means are meant 
to constitute an optimized list. In the next section a short description of the main 
basic concepts involved in any rigorous formulation of the procedure of zeta function 
regularization will be given, together with some results originally obtained by the 
author. 

2. Zeta function of a pseudodifferential operator and determinant 
2.1. The zeta function 

The zeta function £4 of A, a positive-definite elliptic pseudodifferential operator 
('I'DO) of positive order m € R (acting on the space of smooth sections of E, an 
n-dimensional vector bundle over a closed n-dimensional manifold, M) is defined as 

(a(s) = tr A- s = V \J S , Re s > - = s . (4) 

j 

being sq = dimM/ord A the abscissa of convergence of Ca{s). It can be proven that 
Ca{s) has a meromorphic continuation to the whole complex plane C (regular at 
s = 0), provided the principal symbol of A (a m (;r,£)) admits a spectral cut: Lg = 
{A € C; Arg A = 8,8 1 <6 < 6> 2 }, Spec A n Lg = (Agmon-Nirenberg conditiorJUJ) . 
This definition of Ca(s) depends on the position of the cut Lg. The only possible 
singularities of Ca{s) are simple poles at s& = (n — k)/m, k = 0, 1, 2, . . . , n — 1, n + 

1, M. Kontsevich and S. Vishik have managed to extend this definition to the 

case when m G C (no spectral cut exists )P1 
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2.2. Zeta regularized determinant 

Let A be a 'I'DO operator with a spectral decomposition: {ipi, Ai}ig/, with / some set 
of indices. The definition of determinant starts by trying to make sense of the prod- 
uct Yiiei wmcn can be easily transformed into a "sum": lnj} ie/ Xi — lnAj. 
From the definition of the zeta function of A: Ca(s) — J2i£i^7 S > by taking the 
derivative at s = 0: Ca(®) = — we arr i ve to the following definition of 

determinant of A [T2"] : 

det c A = exp[-&(0)]. (5) 

An older definition (due to Weierstrass) is obtained by subtracting in the series 
above (when it is such) the leading behavior of A.; as a function of i, as i — > oo, 
until the series is made to converge [13]. The shortcoming — for physical 

applications — is here that these additional terms turn out to be non-local and, thus, 
are non-admissible in a renormalization procedure. 

In algebraic QFT, to write down an action in operator language one needs a func- 
tional that replaces integration. For the Yang-Mills theory this is the Dixmier trace, 
which is the unique extension of the usual trace to the ideal £( 1 '°°) of the compact 
operators T such that the partial sums of its spectrum diverge logarithmically as the 
number of terms in the sum: o- N (T) = ^^Sq 1 Hj = O(logiV), fiQ > Ui > ■■ ■ The 
definition of the Dixmier trace of T is: Dtr T = lirn/v->oo lo g N &n(T), provided that 
the Cesaro means M(a)(N) of the sequence in are convergent as N — > oo (remem- 
ber that: M(/)(A) = j^j f, f(u)—). Then, the Hardy-Littlewood theorem can be 
stated in a way that connects the Dixmier trace with the residue of the zeta function 
of the operator T _1 at s = 1 (see Connes [13]): Dtr T = lim s _ J . 1 +(s — l)£ T -i(s). 

The Wodzicki (or noncommutative) residuJ^ is the only extension of the 
Dixmier trace to the ^DOs which are not in £( 1,00 ). It is the only trace one can 
define in the algebra of vE'DOs (up to a multiplicative constant), its definition being: 
res A = 2 Res s= o tr (AA~ S ), with A the Laplacian. It satisfies the trace condition: 
res (AB) = res (BA). A very important property is that it can be expressed as an 
integral (local form) res A = J s , M tr a_„(x, £) d£ with S* M C T*M the co-sphere 
bundle on M (some authors put a coefficient in front of the integral: Adler-Manin 
residue^!). 

If dim M = n = — ord A (M compact Riemann, A elliptic, n £ N) it coincides 
with the Dixmier trace, and one has Res s= iC,4 (s) = — res A^ 1 . The Wodzicki residue 
continues to make sense for \&DOs of arbitrary order and, even if the symbols 
aj(x, £), j < to, are not invariant under coordinate choice, their integral is, and 
defines a trace. All residua at poles of the zeta function of a \PDO can be easily 
obtained from the Wodzciki residue^l 

2.3. Multiplicative anomaly 

Given A, B and AB 'J'DOs, even if (a, Cb and (ab exist, it turns out that, in general, 
det^(AB) ^ det^A det^B. The multiplicative (or noncommutative, or determinant) 
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anomaly is defined as: 

det c (AB) 



S(A,B) =ln 



Cab(0)+C'a(0)+Cb(0)- (6) 



_det c A det c B 
Wodzicki's formula for the multiplicative anomaljff^^l 
res{[lna(A,B)] 2 } 

5(A, B) = —± ; i -, a(A, B) := ^ords^-ordA (?) 

v ' ; 2ordAordS(ordA + ordS)' K 1 ' w 

At the level of Quantum Mechanics (QM), where it was originally introduced 
by Feynman, the path-integral approach is just an alternative formulation of the 
theory. In QFT it is much more than this, being in many occasions the actual 
formulation of QFT [T5]. In short, consider the Gaussian functional integration 

j[d$] exp{- J d D x[&{x)( )$(ar) + ...]j — > det ( ) ±1 , (8) 

(the sign ± depends on the spin-class of the integration fields) and assume that the 
operator matrix has the following simple structure (being each A; an operator on 
its own): 

/ a a \ / a \ 

(9) 

where the last expression is the result of diagonalizing the operator matrix. A 
question now arises. What is the determinant of the operator matrix: det(AB) 
or det A ■ det Bl This has been very much on discussion during the last months^. 
There is agreement in that: (i) In a situation where a superselection rule exists, AB 
has no sense (much less its determinant), and then the answer must be det A ■ det B. 
(ii) If the diagonal form is obtained after a change of basis (diagonalization process), 
then the quantity that is preserved by such transformations is the value of det(AB) 
and not the product of the individual determinants (there are counterexamples 
supporting this viewpoinlPH) . 



( A1 


A 2 \ 




Us 


A J 





3. Explicit calculation of Ca and det^ A 

A fundamental property of many zeta functions is the existence of a reflection for- 
mula, also known as the functional equation by mathematicians. For the Riemann 
zeta function: T(s/2)((s) — 7r s ~ 1 / 2 r(l — s/2)£(l — s). For a generic zeta function, 
Z(s), it is Z(u> — s) — F(u), s)Z(s), and allows for its analytic continuation in an 
easy way — what is, as advanced above, the whole story of the zeta function regu- 
larization procedure (at least the main part of it). But the analytically continued 
expression thus obtained is jus t another series, again with a slow convergence be- 
havior, of power series type^ 1 (actually the same that the original series had, in 
its own domain of validity). S. Chowla and A. Selberg found a formula, for the 
Epstein zeta function in the two-dimensional case^l, that yields exponentially quick 
convergence, and not only in the reflected domain. They were extremely proud of 
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that formula — as one can appreciate just reading the original paper (where ac- 
tually no hint about its derivation was given, see [23]). In Ref. [24], I generalized 
this expression to inhomogeneous zeta functions (most important for physical ap- 
plications), but staying always in two dimensions, for this was commonly believed 
to be an unsurmountable restriction of the original formula (see, ej*., Ref. [25]). 
Later I obtained an extension to an arbitrary number of dimensions^, both in the 
homogeneous (quadratic form) and non-homogeneous (quadratic plus affine form) 
cases. 

In short, for the following zeta functions (corresponding to the general quadratic 
— plus affine — case and to the general affine case, in any number of dimensions, d) 
explicit formulas of the CS type were obtained in [26], namely, 



E 

_ rrd 



&(*)= >. [Q(n) + A(n)]- S (10) 

and 



c 2 ( S ) = A (™r s > ( n ) 

aeN d 

where Q is a non-negative quadratic form and A a general affine one, in d dimensions 
(giving rise to Epstein and Barnes zeta functions, respectively). Moreover, expres- 
sions for the more difficult cases when the summation ranges are interchanged, that 
is: 

Ca(«)= E [Q(n)+A(n)}- S (12) 
neN d 

and 

_ rrd 

have been given in [2"5] . 

3.1. Extended Epstein zeta function in p dimensions 

The starting point is Poisson's resummation formula in p dimensions, which 
arises from the distribution identity J2neZ p ^(^ ~ ^) = X)rneZ p e l27rm ' x . (We 
shall indistinctly write rh ■ x = m T x in what follows.) Applying this identity 
to the function f(x) = exp ^— ^x T Ax + b^xj , with A an invertible p x p ma- 
trix, and integrating over x G R p , we get: X^gZ p ex P (j~^n T Ari + b 1 

(27TF 



E™eZ p ex P 



VdetA 

zeta function (Re s > p/2); 



\{b + 2mm\ A~ 1 (b + 27 



Consider now the following 



Ca,«*,00= E 7j(n + c) T A(n + c) + q = ^ ' [Q (n + c) + q] 



aeZ" 
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(the prime on the summation signs mean that the point ft = is to be excluded from 
the sum). The aim is to obtain a formula giving (the analytic continuation of) this 
multidimensional zeta function in terms of an exponentially convergent multiseries 
and which is valid in the whole complex plane, explicitly exhibiting the singularities 
(simple poles) of the meromorphic continuation and the corresponding residua. The 
only condition on the matrix A is that it corresponds to a (non negative) quadratic 
form, which we call Q. The vector c is arbitrary, while q will (for the moment) be 
a positive constant. Use of the Poisson resummation formula yields 

, , (2tt)p/V /2 ~ s T(s-p/2) <2 s / 2 +vl A+2 ii s q- s / 2 +P/ A . . 

CWs) = —jm + — ffir W — (14) 



x cos(27rm • g) {m T A^mf 1 vl K p/2 _ s 2qm T A~^ 

where K v is the modified Bessel function of the second kind, and the subindex 1/2 in 
Z^ 2 means that only half of the vectors m G Z p intervene in the sum. That is, if we 
take an m € Z p we must then exclude —rn (as simple criterion one can, for instance, 
select those vectors in Z p \{0} whose first non-zero component is positive). Eq. (fH| 
fulfills all the requirements demanded before. It is notorious to observe how the only 
pole of this inhomogeneous Epstein zeta function appears explicitly at s = p/2, just 

where it should, its residue being Res s=J ,/2C.A,c,g( s ) — ^/ d ^ tA r(p/2) ' ^ * s re l a tively 
simple to obtain the limit of expression (fH|) as q — > 0. 

When q = there is no way to use the Poisson formula on all p indices of 
ft. However, one can still use it on some of the p indices ft only, say on just 
one of them, n\. Poisson's formula on one index reduces to the celebrated Ja- 
cobi identity for the #3 function, which can be written as X)^L-oo e~^ n+z ^ * = 
^/j 1 + Yl^Li n ^ cos(27rnz) . Here z and t are arbitrary complex numbers, 
z,t E C, with the only restriction that Ret > 0. Applying this last formula to 
the first component, m, we obtain the following recurrent formula (for the sake of 
simplicity we set c = 0, but the result can be easily generalized to c / 0): 

CW«) = C.g,.00 + \fl n S T{s ) /2) Ca^^* - 1/2) + aS/2+ vI r(s) £ I ( 15 ) 

Ef 7Tn itr^\ s -i/2/^t a - , \V4-s/2 f2irn 1 [ZfZ Z 

cos I b n 2 Jn 1 (n 2 A p ^in 2 +q) K s _ 1/2 I y A p _in 2 + q 

ni—l ^ * 

This is a recurrent formula in p, the number of dimensions, the first term of the 
recurrence being (see e.g., the 6th reference in [10j ) 

W-) - 2 E (- 2 + «) s = *~ s + \fl TJf fuf 1 « 1/2 ~ s 

n=l V V ' ' 

+ i^ a - 1 /4- s / 2g i/4- s / 2 £ n*-V*K s _ 1/2 (jhrn^/f) . (16) 
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To take in these expressions the limit q — > is immediate: 



. , s fnT(s- 1/2) , , s 4tt s ' 

C M oW = 2 ^C(2*) + ^- t Ap _ iro>0 (s - 1/2) + a8/2+1/4r(g) £ 

cos I 6 rt 2 ) «i (n 2 A p _iri 2 j a s -i/2 I — -p- y ™2 & P -i n 2 J ■ (17) 

ni- 1 \ V / 



In the above formulas, A is a p x p symmetric matrix A = (ay). ._, 2 = A T , 
Ap^i the (p — 1) x (p — 1) reduced matrix A p ^\ = {a,ij) i ._ 2 , a the component 

a = an, 6 the p — 1 vector b = (021, . . . , a p i) T = (ai2, ■ • ■ , ai P ) T , and A p _i is the 
following (p — 1) x (p — 1) matrix: A p _i = A p _i — j-b®b. 

It turns out that the limit as q — > of Eq. (fT4|) is again the recurrent formula 
(|T7)) . More precisely, what is obtained in the limit is the reflected formula, which 

2s — p/2 

one gets after using Epstein zeta function's reflection T(s)Z(s; A) — 2^==r(p/2 — 

s)Z(p/2 — s; A^ 1 ), being Z(s; A) the Epstein zeta functional This result is easy to 
understand after some thinking. Summing up, we have thus checked that Eq. (|14|) 
is valid for any q > 0, since it contains in a hidden way, for q = 0, the recurrent 
expression (fTT)) . 

The formulas here can be considered as generalizations of the Chowla-Selberg 
series formula. All share the same properties that are so much appreciated by 
number-theoretists as pertaining to the CS formula. In a way, these expressions 
can be viewed as improved reflection formulas for zeta functions; they are in fact 
much better than those in several aspects: while a reflection formula connects one 
region of the complex plane with a complementary region (with some intersection) 
by analytical continuation, the CS formula and the formulas above are valid on the 
whole complex plane, exhibiting the poles of the zeta function and the correspond- 
ing residua explicitly. Even more important, while a reflection formula is intended 
to replace the initial expression of the zeta function — a power series whose con- 
vergence can be extremely slow — by another power series with the same type of 
convergence, it turns out that the expressions here obtained give the meromorphic 
extension of the zeta function, on the whole complex s-plane, in terms of an expo- 
nentially decreasing power series (as was the case with the CS formula, that one 
being its most precious property). Actually, exponential convergence strictly holds 
under the condition that q > 0. However, the formulas themselves are valid for 
q < or even complex. What is not guaranteed for general q eC is the exponential 
convergence of the series. Those analytical continuations in q must be dealt with 
specifically. The physical example of a field theory with a chemical potential falls 
clearly into this class. 
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3.2. Generalized Epstein zeta function in d — 2 



For completeness, let us write down the corresponding series when p = 2 explicitly. 
They are, with q > 0, 



( E {s;a,b,c;q) = -q 3 + 



(!) 



1/4 



(s - i)Va r( s ) 



(18) 



(2tt) s ^V" 1 / 2 cos(7rn6/a) ^V" 2s (A- 



4ag \ 
d 2 / 



l/4-*/2 



A". 



s-1/2 



4aq 



where A = 4ac — 6 2 > 0, and, with (7 = 0, the CS formul 

C £ (s; a, b, c; 0) = 2C(2 S ) a" 5 + f^^m r ( s - 1/2)C(2* - 1) 



r( s )A s - 1 /2 



2 s+5 /2 7rS 



r( s ) A s / 2 -V4 

v 7 v n— 1 



^— = ^ n s - 1 / 2 ( r 1 _ 2s (n) cos(7rn6/a) A s _ 1/2 (™Va) . (19) 



where tr s (n) = X)d|n ^ sum over the s-powers of the divisors of n. We observe that 
the rhs's of (TT8"|) and flT9"|) exhibit a simple pole at s = 1, with common residue: 
Res s= iCB(s;a, 6, c;q) = ^= = Res s=1 ( E (s;a,b,c;0). 



3.3. Truncated Epstein zeta function in d = 2 

The most involved case in the family of Epstein-like zeta functions corresponds 
to having to deal with a truncated range. This comes about when one imposes 
boundary conditions of the usual Dirichlet or Neumann type. Jacobi's theta function 
identity and Poisson's summation formula are then useless and no expression in 
terms of a convergent series for the analytical continuation to values of Re s below 
the abscissa of convergence can be obtained. The method must use then the zeta 
function regularization theorerrJ^ and the best one gets is an asymptotic series. The 
issue of extending the CS formula, or the most general expression we have obtained 
before, to this situation is not an easy one (see, however, Ref. 26 ). This problem 
has seldom (if ever) been properly addressed in the literature. 

As an example, let us consider the following series in one dimension: 
(c(s;a, c;q) = J^L-oo [a( n + c ) 2 + <l] S , Res > 1/2. Associated with this zeta 
functions, but considerably more difficult to treat, is the truncated series, with 
indices running from to oo 

oo 

( Gt (s;a,c;q) = J2H n + c ) 2 +<lT S , Res > 1/2. (20) 

n=0 

In this case the Jacobi identity is of no use. The way to proceed is employing specific 
techniques of analytic continuation of zeta functions. There is no place to describe 
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them here in detail. The usual method involves three steps The first step is 
easy: to write the initial series as a Mellin transform X^^Lo [ a ( n + c ) 2 + 1] — 
TJs) S^Lo Jo°° dtt 3-1 exp { — [a(n + c) 2 + q]t} . The second, to expand in power se- 
ries part of the exponential, leaving a converging factor: J2n°=o [ a ( n + c ) 2 + <?] * = 
rk ^n=o Jo°° dt Eto + c) 2 ^™-^-' 4 . The third, and most difficult, 

step is to interchange the order of the two summations — with the aim to ob- 
tain a series of zeta functions — what means transforming the second series into 
an integral along a path on the complex plane, that has to be closed into a cir- 
cuit (the sum over poles inside reproduces the original series), with a part of it 
being sent to infinity. Usually, after interchanging the first series and the inte- 
gral, there is a contribution of this part of the circuit at infinity, what provides in 
the end an additional contribution to the trivial commutation (given by the zeta 

inoi 

function regularization theorem 1 ^ 1 ). More important, what one obtains in general 
through this process is not a convergent series of zeta functions, but an asymp- 
totic series (see e.g., the 4th and 6th references in [ID]). That is, in our example, 
EZo [«(" + cf + q) s ~ £™ =0 tfrwiffi CiK-am, c) + additional terms. Be- 
ing more precise, as outcome of the whole process we obtain the following result for 
the analytic continuation of the zeta function 



(a. (, «, c; „ ~ - o V ' + £- £ M±fi (i)-<.,-*M (21, 



1 \ _ s q- s (-l) m r(m + s) tq 

2 J Q r(s) ^ m! \a 

-]fl T(S 2 r(l) 2) ql/2 ~ S + ^«- 1/4 - s/2 9 1/4 - s/2 £^- 1/2 cos(2 7 rnc)AVi/ 2 (2 7 rn v ^) 

(Note that this expression re duces to Eq. ([T5]) in the limit c — > 0.) The first se- 
ries on the rhs is asymptotic P^^ . Observe, on the other hand, the singularity 
structure of this zeta function. Apart from the pole at s = 1/2, there is a whole 
sequence of poles at the negative real axis, for s — —1/2, —3/2, . . ., with residua: 
R-es s= i/2_jCG t (s; a, c; q) — j = 0,1,2,... The generalization of this to p 

dimensions can be found in Ref. 1261. 



4. Direct physical applications 

4.1. Calculation of the Casimir energy density 

An application of the pro cedu re is the calc ulation of the Casimir energy density 
corresponding to a masslesJ^or massive^2l scalar field on a general, d dimensional 
toroidal manifold. In the spacetime M = R xE, with S = [0,l] d /^, which is 
topologically equivalent to the d torus, the Casimir energy density for a massive 
scalar field is given directly by Eq. (|14[) at s — —1/2, with q — m 2 (mass of the 
field), 6 = 0, and A being the matrix of the metric g on S, the general c?-torus: 
E%\ m = C g o m 2 ( s = ~ V^)- The components of g are, in fact, the coefficients of the 
different terms of the Laplacian, which is the relevant operator in the Klein-Gordon 
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field equation. The massless case is also obtained, with the same specifications, from 
the corresponding formula Eq. (|17[) . In both cases no extra calculation needs to be 
done, and the physical results follow from a mere identification of the compon ents 
of the matrix A with those of the metric tensor of the manifold in questional 
Very much related with this application but more involved and ambitious is the 
calculation of vacuum energy densities corresponding to spherical configurations 
and the bag model (see [331 131] and references therein). 



4.2. Effective action 



Another application consists in calculating the determinant of a differential opera- 
tor, say the Laplacian on a general p-dimensional torus. A very important p rob lem 
related with this issue is that of the multiplicative anomaly discussed before^. To 
this end the derivative of the zeta function at s = has to be obtained. From Eq. 
dHl) , we get 

j» 4(2aW 4 cos(2vrm-c) „ / , , ^ 



rj (m T A 1 m) 



1/2 



(27r)p/ 2 r(-p/2) g p/ 2 

VdctA 
(-l) k (2Tr) k q k 



p odd, 



[*(fe + l) + 7-lng],p = 2fceven, 

fc! Vdet j4 

and, from here, det A = exp — C^(0). For p — 2, we have explicitly: 



(22) 



— f 00 1 

det A(a, b, c; q) = e M<r-i*<r)/VA (i _ e -**yfiF^ exp i -4 ^ - 



- X! ( 47rnW^ 



cos(irnb/a) d exp 

< / 1 n 



irn / . Aaq 



(23) 



In the homogeneous case (CS formula) we obtain for the determinant: 



det A(a, 6, c) = — exp 



_4 C '(0) - ^ - 4 V co S (7rn6/a)e—^/° 



(24) 



or, in terms of the Teichmiiller coefficients, t% and T2, of the metric tensor (for the 
metric, A, corresponding to the general torus in two dimensions): 



det A(ti,t 2 ) 



T2 

Att 2 \t\ 2 



exp 



-4C'(0) 



7TT 2 
3H^ 



n=l 



V h 



(25) 
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5. Future Perspectives: Operator Regularization 

The Operator Regularization (OR) approach, due originally to D. G. C. McKeon 
and T. N. Sherry^lis considered as a genuine generalization of the zeta regulariza- 
tion approach. Its main aim is to extend zeta regularization, so effective at one-loop 
ordepO-4 t higher loops. It has a distinct advantage over other competing proce- 
dures, in that it can be used with formally non-renormalizable theories, as shown in 
[551 131)] . A further feature of this approach is that divergences are not reabsorbed, 
each one is removed and replaced by an arbitrary factor. Indeed, operator regular- 
ization (OR) does not cure the non-predictability problem of non-renormalizability, 
but an advantage of the method is that the initial Lagrangian does not need to 
be extended with the addition of extra terms. The OR scheme is governed by the 
identity: 



jn 

H- m = lim — 

e^o de n 



1 + (l + aie + a 2 e 2 + ... + a n e n ) — t H 



(26) 



where the a^s are arbitrary, and it is enough that the degree of regularization is 
equal to the loop order, n. 

Two separate aspects of the procedure are, first the regularization itself and, 
second, the analytical continuation, where divergences are replaced by arbitrary 
factors. Thus, the effect of OR is in the end replace the divergent poles by arbitrary 
constants, as 



1 

to yield the finite expression 



(27) 



H m = a n c- n + ■■■+ a x c-x + c . (28) 
5.0.1. Generalization and further extensions 

The OR method can be generalized to multiple operators, as in multi-loop cases 

in 

H - mi . . . H -m r = lim + + a + . . . + «\ 

£^o de n L y ' 

n 

x — H- c - rni ■ ■ ■ (29) 
n\ 

Further extensions of the procedure have been proposed. Let us recall that OR was 
first introduced in the context of the Schwinger approach, 

d n / f™- 1 \ 
InH = - lim — — - H~ e , (30) 
e^o de" V n! / 



which is known to be equivalent to the Feynman one 

>o de n V n\ 



Ml / n 
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The Schwinger form can be transformed into the Feynman one, as 

(m-1)! dH m y ' 

Equivalence with dimensional regularization can be established in many cases, but 
not always. Problems, the main one being unitarity, may appear (see [40]). To start 
with, its naive application to obtain finite amplitudes breaks unitarity. 

A definite advantage of the procedure is that, actually, no symmetry-breaking 
regulating parameter is ever inserted into the initial LagrangiarM^. One can use 
Bogoliubov's recursion formula in order to show how to construct a consistent OR 
operator, and unitarity is upheld by employing a generalized evaluator consistently 
including lower-order quantum corrections to the quantities of interest. Unitarity 
requirements lead to unique expressions for quantum field theoretic quantities, order 
by order in fi. This fact has been proven in many cases (as for the $ 4 theory at 
two-loop order, etc.). But I should say that, to my knowledge, a universal proof of 
this issue is actually still missing. 

A final comment is in order. Using a BPHZ-like scheme, as the above one turns 
out to be, in the end, essentially reintroduces counterterms into the procedure, since 
they are actually hidden in the subtractions taking place at each step. In this way, 
the simplicity of the original zeta function regularization procedure, as described in 
the previous sections, and which is one of its main characteristics, is absent in the 
extended, operator regularization method. 
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